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ABSTRACT 



We have measured the widths of spectral lines from a polar coronal hole using 
the Extreme Ultraviolet Imaging Spectrometer onboard Hinode. Polar coronal 
holes are regions of open magnetic field and the source of the fast solar wind. We 
find that the line widths decrease at relatively low heights. Previous observations 
have attributed such decreases to systematic effects, but we find that such effects 
are too small to explain our results. We conclude that the line narrowing is real. 
The non-thermal line widths are believed to be proportional to the amplitude 
of Alfven waves propagating along these open field lines. Our results suggest 
that Alfven waves are damped at unexpectedly low heights in a polar coronal 
hole. We derive an estimate on the upper limit for the energy dissipated between 
1.1 Rq and 1.3 R & and find that it is enough to account for up to 70% of that 
required to heat the polar coronal hole and accelerate the solar wind. 



1. Introduction 



Despite more than 50 years of investigation, two main aspects of the solar corona 
have not yet been fully explained, namely its heating and the resulting acceleration of the 
solar wind. Two broad classes of models have been proposed for coronal heating and solar 
wind acceleration: wave driven models, which transfer energy to the corona and solar wind 
via waves and turbulence; and reconnection driven models in which energy stored in the 
magnetic fields is releas ed into the cor ona as the fields relax to lower energy states through 



magnetic reconnection ( jCranmer 



20091 ). 



Wave-driven models benefit from the fact that waves have also been detected in coronal 



and interplanetary plasma from the chromosphere to 1 AU (IBanerjee et al. 



20111 ). In 
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particular, Alfven waves appear to be ubiquitous in the Sun. These are transverse MHD 



waves that tra vel along the ma gnetic fiel d lines. They have been detected in chromosp 



bright points (IJess et al. 



solar promi nences ( 



heric 



Okamoto et al. 



solar wind (iBelcher fc Davis 



2009) , spicules (IDe Pontieu et al 



2007 



20071 ). the solar corona (ITomczyk et al 



Mcintosh et a 



20 111 ) 



20071 ). and the 



197ll ). They have also been found to have low frequencies, 



with periods on the order of ~ 100 s, roughly matching the timescales of the observ ed 



photospheric motion s from which they presumably originate (jBahng fc Schwarzschild 



Leighton et al 



1961 



19621 ). And they s eem to carry sufficient 



thereby accelerate the solar wind (IDe Pontieu et al 



20071 ). 



power to heat the corona and 



Strictly speaking, s ome of these observations are 



torsional Alfven waves (IVan Doorsselaere et al. 



probably of fast kink waves , not p ure 



20081 ). However, 



Goossens et al. 



(12009f ) 



emphasizes that kink waves in the corona have mixed properties and a better label would 
be "Alfvenic" . For the discussion of waves throughout this paper we adopt the label Alfven 
to describe the waves, with the caveat that we are using the term broadly. 



Alfven waves are expected to dissipate pri marily 



)y co l 



1991 



isional process es such as 



Cranmer 



20021 ). However, 



viscosity, thermal conductivity, and resistivity (jParkerl 
these processes have characteristic damping lengths of about 2 to 5 R Q . These waves may 
provide the sustained energy deposition required by models to match observed solar wind 
acceleration and proton temperatures. However, predictions indicate that they do not damp 
below ~ 2 R Q and Coulomb collision rates are too low for any energy dissipated above 2 R & 



to be conducted downward to heat the base of the corona. 



Some theoretical work has also shown that Alfven waves can 



i nhomogenous plasma t hrough phase mixing (IHeyvaerts fc Pries 



( jMatthaeus et al. 



19991 ). or resonant absorption (jGoossens et al 



je str ongly damped in an 



1983), turbulent cascade 



201 lh . If such dissipation 



were present, Alfven waves would be a viable mechanism for coronal heating and solar wind 
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acceleration. However, to date there has been no unambiguous observational evidence in 
the solar corona for wave damping at low heights, which we define here as < 1.4 R Q ; and it 
is not clear if the conditions required by these models are actually present in the Sun. 

A common observational approach for detecting the signatures of wave damping is to 
look at the non-thermal broadening of optically thin spectral lines. The unresolved motions 



responsible for this broadening are usuall y attributed to Alfven waves and 



predicted to be 



proportional to the > 


vave a 


2001 


> 


Banerjee et al. 


2009) 



amplitude (e.g., 



Banerjee et al. 



1998 



Doyle et al 



1998 



Moran 



Here, we focus only on observations of polar coronal holes, w hich are main 



Krieger et al. 



y ope n field 



19731 ). One 



regions. Coronal holes are the source region of the fast solar wind (] 
reason for studying these regions is that the fast wind is much less influenced by Coulomb 
collisions than th e slow wind, so that the signatures of other processes at play can be more 



readily observed ( IHollweg 



20081 ) 



Many observations have shown that the line widths in coronal holes initially increase 
with height. T his is also predicted for undamp ed waves based on energy conservation 



considerations ( jHassler et al. 



1990 



Moran 



20011 ) . Some studies, though, have found that 



the widths level off or even decrease abo ye ~ 1.1 R&, w 



damped at heights lower than expected (IBaneriee et al 



2003 



Q'Shea et al 



2005 



Doha fc Solomon 



rich i mplies that the waves are 



1998 



Doyle et al. 



1999 



Moran 



20081 ). Such damping could provide a heating 



source for coronal plasma near where the fast wind begins to be accelerated. However, 
the authors of all these studies have argued that systematic effects have rendered their 
results inconclusive. For example, some of the observations were carried out using the So lar 



Ultra violet Measurement of Emitted Radiation Spectrometer (SUMER; 



Wilhelm et al. 



19951 1 onboard the Solar and Heliospheric Observatory (SOHO). Thes e were s i gnific antly 



affected by instrumental scattered light (described in Section |4.ip and 



Moranl (120031 ) and 
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Dolla fc Solomon! (120081) a t tribut ed the line width decrease to a systematic effect due 



to scattered light. iMoranl ( 120031 ) noticed that lines from ions with different formation 
temperatures behaved differently, which suggests a line-of-sight effect in which the observed 
emission comes from different temperature regions. In another measurement, using 



the Coronal Diagno stic Spectrometer (CDS; 



Harrison et al. 



O'Shea et al. 



19951) also onboard SOHO, 



(120051 ) found a decrease in line widths beginning at about the same height 
where photoexcitation became significant. They suggested there may be a possible 
non-damping physical explanation for the line width decrease, but the exact mechanism 
proposed was unstated. All these potential systematic effects have left unresolved whether 
the observed line width decrease is real or not. 

Here, we present unambiguous evidence that the decrease is real for the line widths in 
polar coronal holes. Our measurements address the uncertainties noted in previous studies. 
We also consider line-of-sight issues and extend the observations to 1.4 R Q , a height larger 
than that reached with either CDS or SUMER. The rest of this paper is organized as 
follows: Section [2] describes the observations and Section [3] describes the analysis methods. 
Our results are reported in Section [4] with a discussion of the various possible systematic 
effects. A possible interpretation of these results in terms of Alfven waves is given in 
Section [5J Section [6] concludes. 



2. Observations 



We have combined data from five observations made with the Extreme Ultraviolet 
Imaging Spectrometer (EIS; 



Culhane et al. 



20071 ) onboard Hinode (jKosugi et al 



20071 ). 



These observations were made on 2009 April 23 at times 12:08, 12:42, 13:16, 13:50, and 
15:17 UT using the 2" slit, which was positioned relative to the central solar meridian 
at, respectively, X = -44.5", -14.5", 15.5", 45.4", and 105.6" (Figured]). In each case 
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the vertical center of the 512" long slit was set at about —1140" so that the height range 
included in the observations extended from about 0.95 R Q to 1.45 R Q . All portions of the 
observations were within the boundaries of the south polar coronal hole. The exposure 
integration time for each pointing was about 30 minutes. 

Images taken the sam e day by the Extreme Ultraviolet Imaging Telescope (EIT; 



Delaboudiniere et al 



19951 ) onboard SOHO were inspected to assess the presence and 



potential importance of plume plasmas in the EIS field of view. The field of view did 
not include any significant plume material at any of the slit positions, and so the EIS 
observations can be considered to consist essentially of only interplume plasma. 

The five observations were averaged together in order to improve the statistical 
accuracy. To do this, each dataset was first processed using the standard EIS analysis 
software to remove spikes, warm pixels, and the CCD dark current. Warm pixels flagged 
by the calibration routine were interpolated using th e recommend ed method, which has 



2010). Systematic drifts in 



been shown to accurately reproduce the line profiles (lYoungi 
the wavelength scale were then corrected. The separate observations were combined by 
averaging over pixels located at the same radius. Finally, the data were binned in the 
vertical direction into bins of 32 pixels each. 



3. Analysis 
3.1. Line Fitting 

We fit Gaussian profiles to the spectrum in order to derive line widths AA from these 
data. Fitting is required to resolve the changes in AA, which here are on the order of 
one tenth of the spectral pixel separation of 0.022 A pixel -1 . We tested the accuracy of 
our fitting procedures by generating synthetic data with a known AA, adding to both the 
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line shape and the background a distribution of random noise, corresponding to that seen 
in the observations, and then running the synthetic line through our fitting procedure. 
Our analysis showed that, on average, the fitting procedure reproduces AA to better than 
about 0.1 mA. However the statistical uncertainty of the individual fit parameters was 
underestimated compared to the standard deviation found by performing the test repeatedly 
with different random errors all drawn from the same distribution. There could be several 
reasons for this. One is that the error bars on the intensity data outside the emission line 
are smaller than the level of background scatter. These intensity data error bars are derived 
from the EIS preparation routines and the weighted averaging from the binning. If these 
errors are too small they could cause some pixels to be weighted too strongly in the fit. 
Another possibility is that the assumed fitting function may not be an exact representation 
of the data. For example, there could be weak lines in the spectrum, while we assume a 
constant or linear background level that does not account for them. 

To account for these issues in the data analysis, we derived line width uncertainties 
from a Monte-Carlo error analysis. After performing an initial fit to the data, normally 
distributed random numbers were added to each data point. These values were scaled so 
that the distribution from which they were generated had a standard deviation equal to the 
residual between each point and the initial fit. Thus, weak features not accounted for by the 
model fitting function are treated as noise. This accounts for systematic errors when the 
fitting function is not a perfect representation of the data. The perturbed data were fit and 
the process was repeated with different random numbers drawn from the same distribution. 
The uncertainties on the parameters were determined by the standard deviation of the 
results from many fits. We found that this gave error bars that were similar in magnitude 
to the scatter seen in the AA results versus radial height R. Thus, this procedure produces 
a more reasonable representation of the actual uncertainties. 
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An additional possible source of uncertainty is that there might be warm pixels that 
are not identified as such in the calibration. The warm pixels maps for the EIS calibration 
are created roughly weekly by identifying those pixels that appear anomalously bright in 
dark frame exposures. These maps show that the warm pixels are distributed randomly 
in both spatial and spectral dimensions. It is possible that there are pixels with a warm 
pixel character that are too weak to appear in the maps. Since any unflagged warm pixels 
are weak by definition it is unlikely that they have a significant effect on the inferred line 
profile. To confirm this, we evaluated their effect on inferred line widths by adding random 
warm pixels to a synthetic spectrum. On average, the spurious pixels tended to slightly 
broaden the lines, but within the fitting uncertainty. In most randomly placed 

warm pixel falls far enough away from the line that it has negligible effect on the inferred 
profile. Because any possible mildly warm pixels fall randomly in space their influence is 
also mitigated by the 32 pixel spatial binning that we used. Thus, the possible existence of 
warm pixels not accounted for in the calibration is not expected to systematically alter the 
inferred line parameters. 



3.2. Line Widths 



The observed width of an optically thin emission line depends on inst rumental 



broad ening; AA; nst , the ion temperature T ; , and the non-thermal velocity v nt ( IPhillips et al 
2008|): 



AA 



FWHM 



AAL t + 41n(2) 



A 



2k B T { 2 



1/2 



(1) 



c J \ M 

Here A is the line wavelength, c is the speed of light, k-Q is the Boltzmann constant, and M 
is the mass of the ion. The thermal plus non-thermal full width at half maximum (FWHM) 
is typically 0.04 - 0.06 A, and the instrumental width is about 0.06 A. 
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The line width includes the sum of the thermal and non-thermal velocities. This 
causes some ambiguity about whether changes in 2] or v nt are responsible for any observed 
radial variation. However, in a polar coronal hole TJ should increase with height due to 
ion cyclotron r esonance heati ng of the ions by waves with frequencies close to the ion 



gyrofrequency (jHollweg 



20081 ). This heating is balanced in the low corona by collisions 



with cooler protons. Previous measurements of 7] are consistent with th e ion temperature 



2009: 


Hahn et al. 


2010 


). At larger h 


rises dramatically 


(Esser et al. 


1999) 



attributed solely to decreases in v nt . Thus, for the purposes of this work it is sufficient to 
consider the effective velocity 



VeS 



2k B Tj 
M 



(2) 



This is the line width expressed as a velocity after subtracting off the instrumental 
broadening. 

The EIS instrumental width is of the same order or even larger than the individual 
thermal and non - therm al line widths. The pre-launch laboratory measurements of 



Korendyke et al. 



( 120061 ) showed that the instrumental width for the 1" slit of EIS was 
0.047 A in the short wavelength b and (170 - 2 1 A) and AA inst = 0.055 A 



in the long wavelength band (240 - 290 A). 



Brown et al. 



( 12008h found that the orbital 



instrumental width is slightly broader with AAi ns t = 0.054 A and 0.055 A in the short 
and long wavelength bands, respectively. A comparison between observations of the same 
location measured with both the 1" and 2" slits showed that the 2" slit has a AAw that is 



0.007 A broader than for the 1" slit ( lYoung 



201 ll ). That is, the instrumental widths in the 



short and long wavelength bands for the 2" slit are 0.061 A and 0.062 A, respectively. 

It is also possible that the instrumental width varies along the length of the slit. This 



has been studied by lYouna (120111 ) and lHara et al.l (120111 ). lYouna (120111 ) measured the line 
widths above the solar limb in equatorial regions and showed that AAi nst has a U-shaped 
dependence on location along the slit. However, because the radial distance from the limb 
varied from rj 1.05 R® at the center of the slit to ~ 1.2 R Q at the ends of the slit, some 
of the apparent variation in AAi ns t may actually have been from radial variation in v nt . 
Additionally, the magnitude of A A w, from this study is larger than expected, being about 



0.066 A at the center of the 2" slit. 



Hara et al. 



( 120 111 ) performed a cross calibration of the 



EIS instrumental width by comparing to ground based observations. They only studied a 



li mited range o 



as 



YouneJ (boilh . 



the EIS detector, but did find qualitiatively the same U-shaped behavior 



For the results presented here we are more concerned with the possible variation of 
AA inst along the slit th a n wit h the absolute value of AA inst , so we use the position- varying 



values given by lYound (|2011| ). which can be accessed through the eis_slit_width routine 
of the EIS analysis software. We have checked the sensitivity of our results to AAi nst 
by deriving v e s assuming the smaller, fixed instrumental widths of 0.061 A and 0.062 A 



discussed above and found that the trends in v e s with height remain essentially the same. 



3.3. Scattered Light: Level 

Superimposed on the true coronal emission is instrumental scattered light, which 
consists of an unshifted ghost spectrum due to scattering of disk radiation by the 
microroughness of the instrumental optics. In our data line widths on the disk were 
narrower than in the off-disk spectra. Thus, scattered light contamination is expected to 
reduce the measured line widths. Far enough above the limb, the instrumental scattered 
light may dominate the measured intensity and cause the observed line width to decrease. 
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There are two issues that must be addressed in order to understand the role of stray 
light in the analysis. First, we must know the actual level of scattered light in EIS. Second, 
we need to determine how the stray light changes the inferred line widths. 

The scattered light contribution to the EI S spectra was tested during an eclipse where 



20101 ) . It was found that the stray 



the moon occulted a portion of the solar disk (lUgartd 
light intensity in the eclipsed portion of the observation was about 2% of the intensity in 
the uneclipsed portion. Previous measurements of the off-disk intensity of lines formed at 



temperatures we 



stray light level ( IHahn et al. 



1 below typical c oronal temperatures have supported this value for the 



201l| ) 



In principle, the amount of stray light in the EIS observations depends on the pointing 
of the instrument as well as on the emission of the solar disk, so that the 2% level found 



by lUgartd (120101 ) may not be representative of the real level of the scattered light in our 
data. To estimate the stray light in our observations we have measured the intensity of a 
chromospheric line. Outside the solar disk, chromospheric lines usually do not emit any 
radiation and the measured intensity is entirely due to scattered light. 

The coldest line observed by EIS that can be used for this estimate is He II 256.23 A. 
This is the only chromospheric line in the EIS spectral range that has enough intensity to 
be observed to large heights off-disk. It is part of a blend with Si x 256.38 A and several 
weaker coronal lines, but the blended lines all lie on the long wavelength side of the He II 
line and can be easily separated using a multi-Gaussian fit. The solid line in Figure [2] shows 
the He II intensity measured in our observation. It is clear that above 1.15 R Q the He II 
intensity is less than 2% of its intensity on-disk. Therefore above 1.15 R & the stray light 
intensity is < 2% of the on-disk intensity. 

Below 1.15 R Q and approaching the solar limb, we observe an increasing level of He II 
intensity. However, this turns out to be real emission, not scattered light. This is because 
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of the large elemental abundance of Helium, which allows there to be a significant He II 
abundance in the corona even though the formation temperature of He II is well below 
the coronal temperature. To estimate the expected coronal intensity of the He II line we 
have calculated the ion charge state dis t ributi on of Helium in an accelerating fast solar 



wind using the model of 



Cranmer et al 



( 120071 ). This model provides the plasma velocity, 



temperature, and density along the entire trajectory of the solar wind starting from the 
lower chromosphere. We have used these values to so l ve the time-dependent equation 



for the charge state distribution following 



Landi et al. 



(j2012al ). We then combined this 



model with the CHIANTI database to calculate the line emissivity at every point along the 
trajectory. Both collisional and radiative excitation were included in the model. For the 
resonant radiative excitation we have used the He II 256.32 A line intensity we measured on 
the solar disk. 

The total line intensity was obtained by integrating the emissivity along the line-of-sight 
assuming a radially expanding coronal hole, and neglecting the pr esence of streamers. W e 



also applied a correction for the difference in density between t he ICranmer et al. 



Hahn et al. 



( 120071 ) 



(120101 ) . This correction reduced 



model and recent density measurements from 
the predicted line intensity by a factor of ~ 1.5. 

The dashed line in Figure [2] indicates the He II 256.32 A intensity predicted by 
the model. Below 1.15 R & the observed emission is almost identical to the predicted 
intensity, both in terms of the absolute intensity an d the rate of decrease. This agreement 
is remarkably good, especially considering that the ICranmer et al.l ( 120071 ) model was 
developed completely independently from our observations and we have used a simplified 
coronal geometry. This excellent agreement implies that even below 1.15 R® the observed 
He II line intensity is not due to scattered light, but rather by real emission ; and that the 



instrumental scattered light is below the 2% threshold estimated by 



IJgartei (120 10f ). Based 



13 



on these results, we conservatively take the stray light level to be 2% of the average on 
disk intensity for the rest of the analysis. We emphasize that this is an upper limit and the 
actual stray light level is likely to be even less. 

3.4. Scattered Light: Effect 

To directly account for the stray light when deriving the line widths, we used a sum of 
two Gaussians to fit the data and extract the line widths. One of the Gaussians in the fit 
had fixed parameters and represented the scattered light profile while the other had free 
parameters and represented the real emission. This procedure is equivalent to subtracting a 
stray light profile from the data. 

First, a single Gaussian fit was performed to the data below the solar limb. The stray 
light line intensity above the limb was set to 2% of the on-disk intensity. The line width 
used for the stray light profile was the same as that extracted from the below limb data, 
but corrected for the varying instrumental width. The centroid position for the stray light 
was the same as that measured on the disk. A two-Gaussian fit to the data above the limb 
was then performed to derive the line width from the real emission. An example of a fit for 
a case with a very high fraction of scattered light is shown in Figure [31 which illustrates 
this procedure for the Fe IX 197.86 A line at about 1.3 R Q . 

In order to estimate the possible influence of scattered light on the inferred line 
widths we also performed measurements using single Gaussian fits to determine the line 
widths and compared them to the two-Gaussian fits described above. The widths from the 
two-Gaussian fits were slightly greater than those from the single Gaussian fits. This is the 
expected result since we essentially subtracted the narrower stray light profiles from the 
data. However, the results from either single or double Gaussian fits agreed to within the 
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uncertainties as long as the stray light fraction was less than ~ 45% of the total intensity. 
For example, for Fe IX at 1.3 R Q , where the scattered light fraction is 42%, the FWHM 
of the Fe IX 197.86 A line was AA = 0.0772 ± 0.0024 A from a single Gaussian fit and 
AA = 0.0803 ± 0.0068 A from the two-Gaussian fit that corrects for stray light (Figure [3]). 

To see the effect of scattered light in more detail we performed single Gaussian fits on 
model spectra. A synthetic spectral line profile was generated as the sum of two Gaussians, 
representing the real and stray light emission. The line widths for the two-Gaussian 
components were set using typical values measured in off-disk and on-disk portions of the 
observation. We also considered small shifts in the centroid position of the scattered light 
component. We then varied the relative intensity of the two components and fit the resulting 
line using a single Gaussian. We found that the output line width varies smoothly between 
the "real" and "stray light" input line widths as the scattered light fraction increases. 
However, in order to produce decreases in the inferred line width of the magnitude found 
in our observations, the scattered light must be > 45% of the total intensity. This confirms 
the result we found from comparing single and double Gaussian fits of the actual data. 

Figure H] shows an example of this model for Fe IX 197.86 A, the line with the largest 
scattered light fraction in the observations. In this example the "real" emission line width 
was set to the maximum observed width, corresponding to 1.12 R & . This overestimates the 
influence of scattered light since the real line width should be increasing with height. Even 
so, the observed line widths decrease more rapidly between 1.1 R Q and 1.3 R Q than can be 
explained by scattered light. 

The line widths described throughout the rest of the paper were determined using the 
two-Gaussian fit to correct for the stray light. In order to ensure that the results we obtain 
are robust to stray light contamination, we adopt a scattered light fraction cutoff of 45% 
and consider only data below this limit. 
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Results 



Effective velocities were derived from the line widths by subtracting the spatially 
varying instrumental width as described in Section 13.21 Figure [5] shows v e s for various 
spectral lines. We focus on these six particular lines because they are strong, relatively 
isolated in the spectrum, and could be measured reliably to large heights. The apparent 



offsets in t uf we attribute to differences in Tj , w 



of the ion ( jLandi fe Cranmer 



2009 



rich is a function of the charge-to-mass ratio 



Hahn et al. 



2010l ). The figure shows that v e s increases 



initially for all ions, but flattens out and begins to decline between 1.1 - 1.2 R@ for Si x, 
Fe IX and Fe x. The Fe XII lines also show some evidence for a decrease, but beginning 
above 1.2 R Q . A clear flattening out is observed for the Fe xm line, though no decrease is 
seen in v e s. 

If the non-thermal velo city is due to Alfven waves then it is re lated to the wave 



amplitude by (5v 2 ) = 2v^ t (jHassler et al. 



1990 



Banerjee et al 



19981 ). The energy flux 



density of the waves is then gi ven by F = p^y^B / ^/ 7 ^ : , wh ere p is the mass density and B 



is the magnetic field strength (IHollweg 



19901 ) 



Moran 



(120011 ) has shown that for undamped 



waves in a flux tube geometry conservation of energy implies that v nt oc p^ 1 / 4 . Thus, in 
a polar coronal hole v nt should increase with R due to the exponential decrease of p with 
height. A constant or decreasing t> nt versus R implies wave damping. 

The dashed lines in Figure [5] roughly illustrate the expected v e s(R) for undamped 
waves based on the p~ 1//4 dependence of v„ t ] where we have e stimated the density from an 



Fe IX intensity ratio and assumed 7] = T e (IHahn et al 



20101 ) . The deviation from these 



trends suggests that all lines do indeed show indications of damping, including the Fe XII 
and Fe xm lines even though their line widths do not decrease as strongly with height. 

Similar results found in the past were mainly ascribed to two effects: instrument- 
scattered light and resonant photoexcitation from radiation emitted by the disk. Here we 
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show that neither of these are an issue with our results. We also consider and rule out 
effects due to the line-of-sight integration implicit in the observations. 



4.1. Scattered Light 

As described in detail in Sections 13.31 and 13. 4[ we accounted for instrument scattered 
light in our observations by performing two-Guassian fits to the data that included a stray 
light profile based on below limb data. We note that when stray light contributes < 45% of 
the total intensity, even the line widths inferred from uncorrected single Gaussian fits were 
found to give the same results. Figure |6] shows the ratio of scattered light to total intensity 
for each of the lines from Figure [5] For each line, stray light makes up less than 10% of the 
total intensity at the heights where the downturn in v e g begins. Thus, even if there were a 
large uncertainty in the stray light intensity, the scattered light fraction would be too small 
to explain the observed behavior in v e ff versus R. 



4.2. Photoexcitation 



O'Shea et al.l (120051 ) showed that the widths of two Mg x lines began to decrease 
at about the same point where photoexcitation became important. They suggested that 
their results were due to resonant photoexcitation systematically affecting the line width 
measurements, although the exact mechanism was unspecified. However, we can rule out 
photoexcitation as an influence on the lines observed here. 

We have measured the relative intensity of the three strong Fe XII lines at 192.39 A, 
193.51 A, and 195.12 A. These are sensitive to resonant scattering of di sk radiation in a 



way similar to the commonly used O VI 1031 A/ 1037 A intensity ratio ( Kohl et al.l 120061 ) . 
These Fe XII lines are due to 3s 2 3p 3A S 3/2 <- 3s 2 3p 2 { 3 P)3d A Pj transitions for J = 1/2, 3/2 



-17- 



5/2, respectively. The Fe XII 193.51 A line blends with an Fe XI line, which prevents it from 
being used in the line width analysis. However, the individual intensities of the two blended 
lines could be determined by estimating the intensity of the Fe XI blend component using 
the measured intensity of the Fe XI 188.30 A line. Both Fe XI lines come from the same 
upper level so their relative intensities are determined solely by a branching ratio. 

Figure [7] shows that the ratios among these Fe XII lines are constant with height. 
If photoexcitation were important the ratios would dec rease with height as the plasma 



19821). Thus, 



changes from being collisionally to radiatively excited (jKohl fe Withbroe 
photoexcitation does not affect the Fe XII lines. 

The other lines in our study do not have such convenient diagnostics for detecting 
photoexcitation, but we can infer that it probably is not a factor for these lines either. First, 
like the Fe XII lines, the Fe x and xm lines start in upper levels connected to the ground 
state by a dipole transition. Compared to the Fe XII lines, the product of the oscillator 



strengths (ILandi et al. 



2012bl ) and the on-disk intensities is about a factor of three less for 
the Fe xm line and about a factor of three greater for the Fe x line. Since these are within 
a factor of a few, the sensitivities of the Fe x and xm lines are probably similar to those 
of the Fe XII lines and we expect that photoexcitation is unimportant. A similar argument 
can be made for the Si x lines used. 

An even stronger argument can be made against photoexcitation influencing the Fe IX 
197.86 A line. This line is due to a 3s 2 3p 5 3d 1 P\ 3s 2 3p 5 Ap 1 So transition. The ground 
state of Fe IX is 3s 2 3p 6 1 S'o. Thus, photoexcitation from the ground state to this upper level 
would involve a strictly forbidden J = - radiative transition. Also, the lower level, 
3s 2 3p 5 3d 1 P i , is connected to the ground level by a strong dipole transition, so that its 
population is very low at any density and radiative pumping from this level to the upper 
1 S'o level is negligible. Since the 3s 2 3p 5 4p 1 So level can essentially only be populated by 
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collisions, photoexcitation is unimportant for this Fe IX line. 



4.3. Line-of-Sight 

Another factor that could influence the observed v e g is the orientation of the magnetic 
field along the line-of-sight of the observations. There are two possible effects. First, 
the outflowing solar wind could have a velocity component along the line-of-sight which 
would produce broadening. Second, the field lines are not always perpendicular to the 
line of sight. Hence fluctuations perpendicular to the magnetic field will be tilted relative 
to the line-of-sight so that the observed v e s appears smaller. Both possibilities affect the 
observations in opposite ways. Because the line intensity is proportional to n^, which falls 
rapidly with height, the data are dominated by the nominal observation height, which is 
the point closest to the Sun and where the line-of-sight is essentially perpendicular to the 
magnetic field. To estimate the magnitude of these effects we assumed that the magnetic 
field lines in the polar coronal hole are radial, that the density scale height is ~ 0.07 R Q 
(corresponding to a typi cal coronal hole temp erature of logT(K) = 6.0) and that the outflow 



velocity is ~ 13 kms 1 (jCranmer et al 



19991 ). With these values we estimate that each of 



these effects may cause line broadening or narrowing on the order of Av e s ~ 0.5 kms -1 . 
This is about 1% the size of the measured values for v e g and thus unimportant. To some 
extent the effects can also be expected to cancel each other out. 



4.4. Differential Emission Measure 

There are some indications that the observation may include multiple plasma 
structures. The profiles for v e s versus R in Figure [5] display a dependence on the ion 
formation temperature. The Fe IX, Fe x and Si x lines form at cooler temperatures than 
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Fe XII or Fe XIII (IBryans et al.ll2009l ) and show an earlier and more pronounced drop in v e g 
with height. To study this in more detail we have performed a Differential Emis sion Measur e 
(DEM ) analysis using the technique described in 



Landi &: Landinil ( 119971 ) and 



Hahn et al 



( 1201 it ). The DEM, <f)(T), shows the distribution of material along the line-of-sight as a 



function of the electron temperature T. 

Figure [8] shows <p(T) at 1.1 R Q . The circles on the plot indicate the points <f)(T t ) 
determined by the measured line intensities, where T t essentially represents the average 
temperature of the plasma from which each emission line originates. The filled circles 
highlight the points corresponding to the lines used to determine v e g(R). The DEM shows 
that the coronal hole emission is dominated by plasma around logT(K) = 6.0, but there is 
a high temperature tail at logT(K) = 6.1 - 6.2. This form of DEM implies that there are 
multiple structures along the line-of-sight. Among the possible interpretations are that the 
cool peak represents emission from the polar coronal hole while the warm tail comes from 
surrounding streamer plasma. Alternatively, the DEM may indicate that there are distinct 
structures with different temperatures within the polar coronal hole. 

The shape of the DEM may explain the observed differences in the v e s(R) behavior of 
the different ions. The fraction of emission originating in the different parts of the DEM 
can be quantified by integrating over the DEM and the contribution function of each line. 
For the purpose of this estimate we take logT(K) = 6.1 to be the dividing line between 
the cool and the warm structure. Integrating from this point up to higher temperatures we 
find that the fraction of emission coming from the high temperature material is 5%, 14%, 
and 43% of the total intensity for Fe IX, Fe x, and Si x, respectively and so these lines 
come primarily from the low temperature peak of the DEM. For Fe XII and Fe xill the 
fraction from the high temperature tail is 67% and 90% of the total intensity, respectively. 
This suggests that the different structures along the line-of-sight could have different v e s(R) 
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profiles as we discuss below. 

Line widths in coronal streamers have been observed to be narrower than in coronal 



holes (IDolla fe Solomon 



20081 ). If the warm material in the DEM represents intervening 



streamer plasma then contamination of the Si x, Fe IX and Fe x by emission from that 
structure could systematically influence the inferred line widths. However, there are 
several factors that indicate that these line widths are dominated by the cool material, 
presumably from the coronal hole. The situation here is analogous to the systematic effect 
of instrumental stray light. In that case we showed that more than 45% of the intensity 
must come from the stray light before the line widths are significantly changed. The 
fraction of light from the warm material is definitely less than 45% for Fe IX and Fe x. 
Additionally, the Si x, Fe IX and Fe x lines all show the same behavior, whereas if there 
were a systematic error from contamination by emission from the warm structure then they 
would behave differently according to the varying level of emission. Therefore, we conclude 
that the velocities derived from Fe IX, Fe x, and Si x reflect the properties of the polar 
coronal hole. 



5. Discussion 



Our analysis shows unambiguous evidence that v nt decreases at relatively low heights 
in a polar coronal hole. These results avoid or resolve the uncertainties that have affected 
previous measurements. Since v nt is thought to be proportional to the Alfven wave 
amplitude, these observations imply that Alfven waves are indeed damped at low heights in 



coronal holes. Previou s studies with S UM 
in v P ff at low heights (IBaneriee et al. 



2005 



Dolla fe Solomo 



n 



1998 



il, CDS, anc 


now 


Dovle et al. 


1999 



)1S have all s 



rown decreases 



Moran 



2003 



O'Shea et al. 



20081 ). Results have been consistent across different instruments 



with varying stray light characteristics and across ions and particular emission lines with 
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Moran 



(120031 ) . where 



varying sensitivity to resonant photoexcitation. This is further evidence that such proposed 
systematic issues are unlikely to explain the observed decrease. It seems that previous 
researchers have overestimated the importance of the systematic un certaint i es on their 
results. The present results are particularly consistent with those of 
an apparent dependence on formation temperature was also noted. 

Based on the DEM analysis, our results imply that Alfven waves are dissipated within 
the cooler structure and that in the warmer structure the damping is less pronounced. If 
the warmer structures are interpreted as streamer plasma along the line-of-sight, our results 
indicate that Alfven wave damping is active at low heights in the coronal hole and hence 
in the fast wind, but is not as strong in streamer material where the slow wind is believed 
to originate. If Alfven wave damping occurs via turbulent cascades of wave power from 
low frequencies to frequencies high enough to cause ion-cyclotron acceleration, our results 
are consistent with signa tures of ion-cyclotron effects being measured in coronal holes 



( ILandi fc Cranmer 



20091 ) and not in streamers (ILandi 



20071 ) 



Though dissipation of Alfven waves by vis cosity, thermal conductivity, and resistivity 



1991 



Cranmer 



20021 ) . damping may 



are not expected to be important below 2 R Q (jParkerl 
be enhanced in an inhomogenous plasma. For example, it has been shown that phase-mixing 
induced by adjacent Alfven waves on neighboring field lines resonat ing out of phase wi th 



each other can enhance the viscous and resistive dissipation rates (jHeyvaerts & Priest 



19831 ). Or the interaction of outward and inward propagating Alfven waves can lead to 



a turbulent cascade that transfers energy to hig her frequencies and for ms other types of 



plasma waves, which are more strongly damped (IMatthaeus et al. 



19991 ). 



Regardless of the physical mechanism that leads to the wave dissipation, the energy 
must be deposited in the low corona and contribute to coronal heating. We can estimate 
an upper bound for the dissipated energy using the expression for the energy flux of an 
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Alfven wave F = 2p v^Va (IDoyle et al. 



19981 ). where Va = B/y/Anp is the Alfven speed. 
The electron density was measured from an Fe IX line intensity ratio to be ~ 6 x 10 7 cm -3 
at 1.1 Rq. Due to the low intensity, the ratio could not be used directly above ~ 1.15 R Q , 
but extrapolating with a scale- height exponential decay gives n B ~9x 10 6 cm" 3 at 1.3 R Q . 



1998); and 



This value agrees with measurements for polar coronal holes ( IWilhelm et al 
since F varies as y/n^, the uncertainty incurred by extrapolating the density should not 
seriously affect our estimate. To estimate the spatially varying magnetic fi eld strength for 



a sup er-radially expanding polar corona 



1999|) and take B w 7 G at 1 R Q flWang 



corona 


hole 


(Wang 


2010) 



hole we used an empirical model ( jCranmer et al. 



Va(1-3-R q ) ~ 1600 kms l . We obtain a n upper bound on 



the ion temperature of 7] ~ T e pa 10 6 K (ILandi fe Cranmer 



v„t u s ing the lower bou nd on 



2009 



Hahn et al. 



2010|). Using 



Equation 2 we find from the Fe x lines that v nt = 42.6 kms -1 and 41.8 kms" 1 at 1.1 Rq 
and 1.3 R & , respectively. Similar results are seen for Fe IX and Si x. Thus, we obtain upper 
bounds of F — 5.4 x 10 5 erg cm -2 s" 1 at 1.1 R Q and F = 9.2 x 10 4 erg cm -2 s" 1 at 1.3 R Q . 
As described in Section EJ there is some uncertainty in the EIS instrumental line width and 
using the smaller estimate would change these values to F = 7.0 x 10 5 erg cm -2 s" 1 and 
F — 1.1 x 10 5 erg cm" 2 s" 1 at 1.1 R Q and 1.3 R & , respectively. This represents an energy 
flux loss of 5.9 x 10 5 erg cm" 2 s . The amount of energy required to he at the coronal hole 



and a ccelerate the fast solar wind is estimated to be 8 x 10 erg cm s ( IWithbroe h Noyes 



19771 ). Thus, we estimate that in our observations Alfven waves may deposit up to about 



70% of the energy required to heat the coronal hole and accelerate the fast solar wind. The 
challenge for the field is now to derive an observational lower limit for the energy deposited 
by Alfven waves in coronal holes. Only then will we be able to determine the fraction of 
heating due to Alfven waves and that due to magnetic reconnection. 
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6. Summary 

We have measured the variation of spectral line widths from 1.05 - 1.40 R Q over a 
polar coronal hole. These line widths are thought to be proportional to the Alfven wave 
amplitude. We found that they deviate from the predicted p _1//4 dependence for undamped 
waves. We have investigated possible systematic effects such as instrumental scattered 
light, resonant photoexcitation, and line-of-sight observational effects and determined that 
they are all too small to explain our observations. Thus, our results suggest that Alfven 
waves are damped at surprisingly low heights in a polar coronal hole. We estimate that the 
amount of energy dissipated can account for a large fraction of that required to heat the 
coronal hole and accelerate the solar wind. 

We thank the anonymous referee for helpful suggestions. MH and DWS were supported 
in part by the NASA Solar Heliospheric Physics program grant NNX09AB25G and the NSF 
Division of Atmospheric and Geospace Sciences SHINE program grant AGS-1060194. The 
work of EL is supported by the NASA grants NNX10AQ58G and NNX11AC20G. 



Table 1:: Line List for DEM Analysis. 



Ion 


A (A) 1 




Transition 1 


Mg VI 


268.991 


2s 2 2p 3 2 D 3/2 


- 2s2p 42 P 1/2 


Mg vi < 


C 270.391 


2s 2 2p 32 D 5/2 


2s2p i2 P 3/2 




L 270.400 


2s 2 2p 3 2 D 3/2 


2s2p 42 P 3/2 


Mg vii 


276.154 


2s 2 2p 2 3 P 


- 2s2p 33 S 1 


Si vi 


249.125 


2s2 2p5 2 P 1/2 


- 2s 2pQ 2 Si /2 


Si vii 


272.648 


2s 2 2p i3 P 2 


- 2s2p 53 P l 


Si vii 


275.361 


2s 2 2p 4 3 P 2 


- 2s 2p 5 3 P 2 
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Si vii 


275.676 


2s 2 2p 53 Pi 




2s2p 53 Pi 


Si ix 


258.082 


2s 2 2p 2 1 D 2 




2s 2p 3 1 P> 2 


Si x 


258.374 


2s 2 2p 2 P 3/2 




2s 2p 2 2 P,/ 2 

^ x o/z 


Si x 


261.057 


2s 2 2» 2 P 3/2 




2s 2» 2 2 P /2 


Si x 


271.992 


2s 2 2n 2 Pi/9 

1 1/z 




2s2» 22 S'i/2 


Si x 


277.264 


2s 2 2» 2 P, /9 




2s2» 22 S'i/2 


Fe viii 


185.213 


3» 6 3d 2 Dcjo 




3» 5 3d 2 f 3 P) 2 P 7 /o 


Fe viii 


186.599 


3v 6 3d 2 P> 3/2 

Sr ~^o/Z 




3p 5 3d 2 ( 3 P) 2 P 5/2 

^'/^ V / O/Z 


Fe viii 


194.661 


3» 6 3rf 2 P 5 / 2 




3p 6 4 p 2p 

o/Z 


Fe IX 


188.497 


3s 2 3p 5 3rf 3 P 4 




3s 2 3p 4 ( 3 P) 3rf 23 G 5 


Fe IX 


189.941 


3s 2 3p 5 3rf 3 P 3 




3s 2 3p 4 ( 3 P) 3rf 23 G 4 


Fe IX 


197.862 


3s 2 3p 5 S^Pi 


— 


3s 2 3p 5 4p 1 1 S 


Fe x 


174.531 


3s 2 3» 5 2 Pi /2 




3s 2 3p 4 ( 3 P) 3d 2 D 5/2 


Fe x 


184.537 


3s 2 3» 52 P^/2 

7 «-* J t o/z 




3s 2 3p 4 ( X P) 3d 2 Si /2 


Fe x 


190.037 


3s 2 3p 5 2 P io 

J u -*- 1/z 




3s 2 3p 4 ( X P) 3d 2 Si/ 2 


Fe x 


193.715 


3s 2 3n 52 P./2 

7 ' - 1 o/z 




3s 2 3p 4 ( 1 S') 3d 2 D 5/2 


Fe x 


t 257.259 


3s 2 3n 52 P,/2 

J *-* J l' - 1 o/z 




3s 2 3p 4 ( 3 P) 3d 4 D 5/2 








^ 257.263 


3s 2 3p 52 P^/2 

7 - 1 o/z 




3s 2 3p 4 ( 3 P) 3d A D 7/2 


Fe XI 


180.401 


3s 2 3p 43 P 2 




3s 2 3p 3 ( 4 S) 3d 3 P> 3 


T7V-\ "v t 

re -X.1 


1 SO 1 R7 
loz,. 10 ( 


OS 0J9 J i 




q„2 q„3 /4 m Qj3n 
os op ^ O J 00 ±^2 


Fe XI 


188.217 


3s 2 3p A 3 P 2 




3s 2 3p 3 ( 2 D)3d 3 P 2 


Fe XI 


188.299 


3s 2 3p 4 3 P 2 




3s 2 3p 3 { 2 D) 3d x Pi 


Fe XI 


189.711 


3s 2 3p 43 P 




3s 2 3p 3 ( 2 P)3rf 3 Pi 


Fe XII 


192.394 


3s 2 Sp 34 ^ 




3s 2 3p 2 ( 3 P) 3d 4 P 1/2 


Fe XII 


193.509 


3s 2 Sp 34 ^ 




3s 2 3p 2 ( 3 P) 3d 4 P 3/2 
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* Fe XII 


195.119 


3s 2 3p 34 S 3/2 


- 3s 2 3p 2 ( 3 P)3d 4 P b/2 


* Fe xm 


202.044 


3s 2 3p 2 3 P 


- 3s 2 3p3d 3 P 1 




j- 203.797 


3s 2 3p 2 3 P 2 


- 3s 2 3p3d 3 D 2 


Fe xm < 






203.827 


3s 2 3p 2 3 P 2 


- 3s 2 3p3d 3 D 3 


Fe xm 


251.953 


3s 2 3p 2 3 P 2 


- 3s 3p 3 3 S 1 


Fe xiv 


264.789 


3s 2 3p 2 P 3/2 


3s3p 22 P 3/2 


Fe xiv 


270.521 


3s 2 3p 2 P 3/2 


- 3s3p 22 P l/2 


Fe xv 


284.163 


3s 2 % 


- 3s 3p x Pi 



1 Wavelengths and transitions taken from CHIANTI (IDere et al.lll997l : lLandi et al.ll2012bl ). 
°Brackets indicate blends from the same ion. 
Asterisks mark the lines used in the analysis of the line widths. 
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Fig. 1. — The slit positions of the EIS observations are illustrated here overlayed on an 
EIT/ SOHO image in the 195 A band, which consists primarily of Fe XII emission. 
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Fig. 2. — The solid line shows the measured He II 256.32 A line intensity in the observation. 
The horizontal dotted line indicates indicates the intensity that is 2% of the on-disk average. 
Above about 1.15 R Q the measured intensity is clearly less than this value, which is therefore 
an upper bound on the scattered light. The dashed line illustrates the model prediction for 
the real emission from the He II line. We expect the difference between the measurement 
and the model to be caused by instrumental scattered light. The close agreement between 
the measured and predicted emission implies that the scattered light level is below 2% for 
all of the off-limb data. 
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197.70 197.80 197.90 198.00 

Wavelength (A) 



Fig. 3. — Two-Gaussian fit to the Fe IX 197.86 A line at 1.3 R®. Here the dotted line is the 
fixed stray light profile with a FWHM 0.074 A, the dashed line is the inferred profile from 
real emission giving FWHM 0.0803 ± 0.0068 A, and the solid line is the sum of the two. A 
single Gaussian fit to the same data would give AA = 0.0772 ± 0.0024 A. This is an extreme 
example with a very high stray light level of 42%. 
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Fig. 4. — The thick solid line shows the Gaussian line width that would be inferred from a 
single Gaussian fit as a function of the scattered light fraction. In this model the line width 
of the real component is set to the width of the Fe IX 197.86 A line at its maximum at 
1.12 R Q . The line width of the scattered component is set to the width of the same line on 
the solar disk. The points on the plot show the observed line width of the Fe IX line with 
labels on selected points to indicate the corresponding height in units of R Q . The dashed 
line is a linear fit to the observations between 1.1 R Q and 1.3 Rq to illustrate that the lines 
narrow more rapidly than predicted by the model when scattered light is less than ~ 45% of 
the total line intensity. The la uncertainty in the slope of the fit is indicated by the dotted 
lines. 
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1.0 1.1 1.2 1.3 1.4 

Distance from Sun Center (R Q ) 



Fig. 5. — Line widths, expressed as f e g versus height. Points where the instrumental scat- 
tered light is greater than 45% have been omitted. The dashed lines illustrate, for each line, 
the f e ff expected for undamped waves normalized at 1.1 R e . 
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• Si X 261.04A 




1.0 1.1 1.2 1.3 1.4 

Distance from Sun Center (R Q ) 



Fig. 6. — Scattered light intensity in each line as a fraction of the total intensity for the line 
widths shown in Figure [5j The scattered light fraction for each line is below about 10% at 
the heights where damping begins to be seen in Figure [51 
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Fig. 7.— The intensity ratios of Fe xn 3s 2 3p 3i S 3/2 <- 3s 2 3p 2 ( 3 P)3d 4 Pj for J =1/2, 3/2, 
and 5/2, corresponding to the 192.39 A, 193.51 A, and 195.12 A lines, respectively. These 
line ratios are sensitive to resonant photoexcitation. The dashed lines show the average ratio 
between 1.05 and 1.15 R Q . If photoexcitation we re significant these ratio s would decrease 



with height as the collisional excitation rate drops (IKohl fe Withbroe 



19821 ). That the ratios 



are constant with height demonstrates that photoexcitation does not affect these Fe XII lines. 
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Fig. 8. — The DEM at 1.1 R Q . The open and filled circles show the points based on the 
measured line intensities used to calculate the DEM. The scatter in these points gives an 
estimate of the uncertainty in the DEM. The filled circles highlight those particular points 
corresponding to the lines in Figure [5] 
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